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•.;/- SYSTEM AND METHOD FOR AIR FLOW AND EGR FLOW ESTIMATION 

5 Cross Reference to Related Applications 

This application is a divisional of Serial No. 09/539,227 
filed March 30, 2000, which is a continuation-in-part of 
Application Serial No. 09/226,681, filed January 11, 1999 and 
Application Serial No. 09/517,201, filed March 2, 2000, which is, 

10 in turn, a continuation-in-part of Serial No. 09/226,681, filed 
January 11, 1999, and Serial No. 09/090,125, filed June 4, 1998. 
This application hereby expressly incorporates by reference the 
entire specification of Serial No. 09/226,681, filed 
January 11, 1999, Serial No. 09/517,201 filed March 2, 2000, and 

15 Serial No. 09/090,125, filed June 4, 1998. 

Field of the Invention 
The present invention relates to an air/fuel ratio control 
system for an internal combustion engine where air flow and 
2 0 exhaust gas recirculation flow are calculated from pressure 
sensors . 

Background of the Invention 
Engine control systems often determine the amount of fuel to 
25 inject by measuring a manifold pressure, along with other engine 
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operating conditions. This method is often referred to by those 
skilled in the art as the speed density method. In this method, 
a mean value model of engine operation is constructed, where an 
? .. average manifold pressure at a given speed results in a certain 

5 air flow into the cylinder. In this type of system, measurement 
of the manifold pressure is critical for proper prediction of the 
air flow into the cylinder and thus for proper air/fuel ratio 
control . 

As stated above, many methods are available to estimate 
10 cylinder air charge using a manifold pressure sensor. Typically, 
engine maps are provided that provide a cylinder air charge as a 
•y,\ function of measured manifold pressure, manifold temperature, and 

engine speed. In engines that also utilize exhaust gas 

: -t. 

[..,y. recirculation, an improved cylinder air charge estimate is 

:;^-:£l'5 obtained by providing adjustments based on the amount of exhaust 
• : gas recirculation. 

One particular method is described in U.S. Patent 5,205,260. 
In this method, an EGR flow is estimated based on differential 
pressure across a flow control valve and based on a cross- 
20 sectional area of the valve. Then, this flow is used in a 

manifold filling model to estimate a partial pressure of EGR in 
the intake manifold. Then, based on this partial pressure of EGR 
arid measured manifold pressure, a cylinder air charge value is 
computed. 
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The inventors herein have recognized a disadvantage with the 
above system. In particular, estimating EGR flow in this manner 
leads to estimation inaccuracies, due to valve area uncertainty. 
This uncertainty may be caused by deposits in the valve. Since 

5 inaccuracies in EGR flow directly affect estimated cylinder air 
charge, this leads to inaccuracies in calculated fuel injection 
amount and therefore may degrade air- fuel ratio control. 

Another approach to determining EGR amount has used a 
differential pressure measurement across an orifice to infer a 

0 flow of exhaust gas. Traditionally, the orifice is located 
upstream of the exhaust gas recirculation flow control valve. 
Thus, the pressure measurements are shielded from the intake 
manifold pressure pulsations; however, the pressure measurements 
are not shielded from the exhaust pressure pulsations. In the 

5 traditional system, the high frequency pressure pulsations 
present in the pressure measurements are reduced by using a 
conventional low pass filter. Such a system is disclosed in U.S. 
.Patent 5,613,479. 

The inventors herein have recognized a significant 

0 opportunity to reduce total system cost by relocating the orifice 
downstream of the exhaust gas recirculation flow control valve 
but before the intake manifold. Thus, the manifold pressure 
sensor can be used to measure the pressure downstream of the 
orifice and a single absolute pressure sensor can be used to 

5 measure the pressure upstream of the orifice. This creates the 
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needed differential pressure to measure exhaust gas recirculation 
flow, as well as an opportunity to estimate cylinder fresh 
charge . 

5 Summary of the Invention 

An object of the invention claimed herein is to provide an 
exhaust gas recirculation measurement system and cylinder air 
charge estimation system with improve accuracy. 

The above object is achieved, and problems of prior 
•:10 approaches overcome, by a system for estimating engine flows, 
including exhaust gas flow from an exhaust manifold of an 
internal combustion engine to an intake manifold of the engine, '■ 
the system comprising: a flow control valve having a variable 
orifice positioned in an exhaust gas recirculation path between 
.15 the exhaust manifold and intake manifold of the engine; a fixed 
orifice area located in said path and downstream of said valve; 
and a computer for measuring a first pressure between said fixed 
orifice area and said flow control valve, measuring a second 
' pressure downstream of said fixed orifice area, calculating a 
2 0 recirculated exhaust flow based on said first pressure and said 
second pressure, and determining an air amount based on said 
second pressure and said calculated recirculated exhaust flow. 

By using common signals for estimating EGR flow and cylinder 
air charge, a simplified structure and reduced cost system is 
25 obtained. Further, by using a fixed orifice area downstream of a 



- 4 - 



control valve to recirculate exhaust gas along with pressure 
measurements upstream and downstream of the orifice, a more 
accurate EGR flow estimate is obtained. In particular, the 
estimation scheme does not have to account for the changing valve 
5 area, and thus less affects have to be included. Further yet, 
the present invention does not need to measure or infer exhaust 
manifold temperature or exhaust manifold pressure. 

An advantage of the above aspect of the invention is that 
more accurate cylinder air charge estimate is obtained. 
10 Another advantage of the above aspect of the invention is 

that the more accurate cylinder air estimate yields improved 
emissions at a reduced system cost. 

Other objects, features and advantages of the present 
■invention will be readily appreciated by the reader of this 
15 specification . 

Brief Description of the Drawings 
The object and advantages described herein will be more 
fully understood by reading an example of an embodiment in which 
2 0 the invention is used to advantage, referred to herein as the 
Description of the Preferred Embodiment, with reference to the 
drawings wherein: 

Figure 1 is a block diagram of an engine in which the 
invention is used; 
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Figures 2A and 2B are alternate embodiments of the present 
invention; 

Figures 3-6 are a high level flowcharts of various routines 
for controlling EGR flow; 
5 Figure 7 is a block diagram of an engine in which the 

invention is used; 

Figures 8-10 are high level flowcharts of various operations 
performed by a portion of the embodiment shown in Figure 8; 

Figures 11-12 are examples of a fluctuating waveform on 
10'-. which the invention is used to advantage; and 

Figures 13A and 13B are plots showing frequency content of a 
pressure signal and an example of a notch filter's magnitude 
frequency characteristics. 

15 Description of an Embodiment 

Internal combustion engine 10 comprising a plurality of 
cylinders, one cylinder of which is shown in Figure 1, is 
controlled by electronic engine controller 12. Engine 10 
includes combustion chamber 30 and cylinder walls 32 with piston 

20 36 positioned therein and connected to crankshaft 40. Combustion 
". chamber 3 0 communicates with intake manifold 44 and exhaust 

manifold 48 via respective intake valve 52 and exhaust valve 54. 
Exhaust gas oxygen sensor 16 is coupled to exhaust manifold 48 of 
engine 10 upstream of catalytic converter 20. 
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Intake manifold 44 communicates with throttle body 64 via 
throttle plate 66. Intake manifold 44 is also shown having fuel 
injector 68 coupled thereto for delivering fuel in proportion to 
the pulse width of signal (fpw) from controller 12. Fuel is 
5 delivered to fuel injector 68 by a conventional fuel system (not 
shown) including a fuel tank, fuel pump, and fuel rail (not 
shown) . Engine 10 further includes conventional distributorless 
ignition system 88 to provide ignition spark to combustion 
chamber 30 via spark plug 92 in response to controller 12. In 

10' the embodiment described herein, controller 12 is a conventional 
microcomputer including: microprocessor unit 102, input /output 
ports 104, electronic memory chip 106, which is an electronically 
programmable memory in this particular example, random access 
memory 108, and a conventional data bus. 

15 Controller 12 receives various signals from sensors coupled 

to engine 10, in addition to those signals previously discussed, 
including: measurements of inducted mass air flow (MAF) from mass 
■ air flow sensor 110 coupled to throttle body 64; engine coolant 
temperature (ECT) from temperature sensor 112 coupled to cooling 

20 jacket 114; a measurement of manifold pressure (MAP) from 

manifold pressure sensor 116 coupled to intake manifold 44; a 
measurement of throttle position (TP) from throttle position 
sensor 117 coupled to throttle plate 66; and a profile ignition 
pickup signal (PIP) from Hall effect sensor 118 coupled to 

25 . crankshaft 40 indicating and engine speed (N) . 
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Exhaust gas is delivered to intake manifold 44 by a 
conventional EGR tube 202 communicating with exhaust manifold 48, 
EGR valve assembly 200, and EGR orifice 205. Alternatively, tube 
202 could be a internally routed passage in the engine that 
5 communicates between exhaust manifold 48 and intake manifold 44. 
Flow Sensor 2 06 communicates with EGR tube 2 02 between valve 
assembly 200 and orifice 205. Flow sensor 206 also communicates 
with intake manifold 44. Stated another way, exhaust gas travels 
from exhaust manifold 44 first through valve assembly 2 00, then 
10 through EGR orifice 205, to intake manifold 44. EGR valve 

assembly 200 can then be said to be located upstream of orifice 
205. 

Flow sensor 2 06 provides a measurement of manifold pressure 
(MAP) and pressure drop across orifice 205 (DP) to controller 12. 
15 Signals MAP and DP are then used to calculated EGR flow as 

described later herein with particular reference to Figure 3-5. 
EGR valve assembly 2 00 has a valve position (not shown) for 
. controlling a variable area restriction in EGR tube 202, which 
thereby controls EGR flow. EGR valve assembly 200 can either 
20 minimally restrict EGR flow through tube 202 or completely 

restrict EGR flow through tube 202. Vacuum regulator 224 is 
coupled to EGR valve assembly 200. Vacuum regulator 224 receives 
actuation signal (226) from controller 12 for controlling valve 
- position of EGR valve assembly 200. In a preferred embodiment, 
25 EGR valve assembly 200 is a vacuum actuated valve. However, as 
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is obvious to those skilled in the art, any type of flow control 
valve may be used, such as, for example, an electrical solenoid 
powered valve or a stepper motor powered valve. 

Referring now to Figures 2A and 2B, and in particular to 
5 Figure 2A, an alternative embodiment of the present invention is 
shown in which housing 250 contains path 252 with inlet end 254 
and outlet end 256. Variable orifice 258 is controlled by pintle 
260 of valve 200. Housing 250 also holds vacuum regulator 224 
which is coupled to valve 200 and thereby regulates pintle 260. 

10 Path 252 also has orifice 205" coupled to outlet end 256. 

Differential pressure sensor 262 measures pressure difference 
across orifice 205 and provides differential pressure signal 266 
to circuit 268. Pressure sensor 264 measures communicates via 
measurement path 269 with outlet end 256 and measure pressure 

15 downstream of orifice 205 and provides pressure signal 270 to 
circuit 268. Circuit 268 calculates, either digitally using 
microprocessor circuits known to those skilled in the art or 
using analog circuits known to those skilled in the art, the 
product of signals 266 and 270. Circuit 268 then makes the 

20 result of this calculation available in signal 272. 

Alternatively, as shown in Figure 2B, differential sensor 
262 and sensor 264 communicate with downstream flow (not shown) 
via second communication path 274. In this embodiment, paths 256 
and 274 are adapted to be connected to an intake manifold of an 

25 internal combustion. Then, path 274 and 256 will be in fluid 
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communication via the intake manifold. Such an arrangement is 
preferable if circuit 268 also provide signal 276 representing 
the pressure measured by sensor 264. 

Referring now to Figure 3, a routine for calculating EGR 
5 flow (EM) is described. In step 210, the signal MAP is read by 
controller 12 from sensor 2 06, giving a measure of pressure 
downstream of orifice 205. Then, in step 212, the differential 
pressure, DP, across orifice 205 is read by controller 12 from 
sensor 206. In step 214, a correction factor, CF1, partially 

10 accounting for the compressibility effects of the EGR flow is 
calculated as the absolute pressure measured by signal MAP. 
Alternatively, if the downstream pressure measured in step 210 
was pressure relative to atmosphere, correction factor CF1 would 
be calculated as the sum of the pressure relative to atmosphere 

15 plus the absolute pressure due to the atmosphere. Then, in step 
216, EGR flow, EM, is calculated as the square root of the 
product of correction factor CF1, differential pressure DP, and 
constant K. Constant K represents a calibration term that 
accounts for various unit conversions and the area of orifice 

20 205. In this way, pressure and temperature effects due to the 
expansion of the EGR flow through valve 200 are sufficiently 
removed and measurement error is reduced. 

The routine described in Figure 3 exploits the nature of the 
flow due to expansion first through flow control valve 2 00 and 

2 5 then through orifice 2 05, where the source of flow is exhaust 
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manifold 48 and the sink is intake manifold 44 of internal 
combustion engine 10. Due to the typical ranges of exhaust 
manifold pressure and temperature and intake manifold pressure 
(MAP) , EGR flow may be approximated using the product of pressure 
5 difference (DP) across orifice 205 and pressure downstream (MAP) 
of orifice 205 without need for measuring temperature upstream of 
orifice 205 (downstream of flow control valve 200) . 

Referring now to Figure 4, an alternate routine for 
calculating EGR flow (EM) is described. In step 310, the signal 

10 MAP is read by controller 12 from sensor 206, giving a measure of 
pressure downstream of orifice 205. Then, in step 312, the 
differential pressure, DP, across orifice 205 is read by 
controller 12 from sensor 206. In step 314, a correction factor, 
CF1, partially accounting for the compressibility effects of the 

15 EGR flow is calculated as the absolute pressure measured by 

signal MAP. Alternatively, if the downstream pressure measured 
in step 310 was pressure relative to atmosphere, correction 
factor CF1 would be calculated as the sum of the pressure 
relative to atmosphere plus the absolute pressure due to the 

2 0 atmosphere. Then, in step 316, correction factor CF2 is 

calculated as a function of both differential pressure DP and 
downstream pressure MAP, where k represents the ratio of specific 
heats of exhaust gas. Correction factor CF2 further accounts for 
the compressibility effects of the EGR flow. Then, in step 318, 

25 correction factor CF3 is calculated as a function of flow through 
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the engine, MAF . Correction factor CF3 accounts for variations 
in exhaust pressure. Function h represents a function relating 
airflow through the engine (MAF) to exhaust pressure and is 
determined experimentally. Additionally, function h can include 
5 a correction for barometric pressure. In other words, the 
.exhaust pressure is calculated as a function of both MAF and 
barometric pressure. The effect of barometric pressure on 
exhaust pressure is also determined experimentally. Barometric 
pressure can be either measured or estimated using methods known 

10 to those skilled in the art. Then, in step 320, EGR flow, EM, is 
calculated as a function of correction factors CF1, CF2 , CF3, 
differential pressure DP and constant K. In this way, pressure 
and temperature effects due to the expansion of the EGR flow 
through valve 200 are further removed and measurement error is 

15 . further reduced with additional complexity. 

Referring now to Figure 5, a routine for controlling EGR 
flow is described. In step 410, the desired EGR flow, DESEM, is 
calculated as a function of engine operating conditions that 
include engine speed (determined from signal PIP) and airflow 

20 (MAF) . Then, the value of EM calculated according to either 

Figure 3 or 4 is subtracted from DESEM to create an error signal, 
ERROR. Then, in step 414, actuation signal 226 is calculated as 
a function (f) of signal ERROR. In a preferred embodiment, 
function (f) represents a PID controller. Alternatively, 
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function (f) may represent any type of feedback or feedforward 
controller known to those skilled in the art. 

Referring now to Figure 6, a routine for calculating 
cylinder air charge. The routine is executed once per engine 
5 firing. In other words, the routine is executed synchronously 
with engine firing events. First, in step 610, total cylinder 
mass charge (which is the sum of cylinder air amount and cylinder 
exhaust gas recirculation amount) is determined based on manifold 
pressure and manifold temperature. In a preferred embodiment, 
10 the following equation is used where slope (s) and offset (o) are 
determined as a function of engine speed. 
MAP * s~o 

™cyl = j 

m 

Next, in step 612, EGR flow EM divided by engine speed (N) and 
number of cylinders (n cyl ) is determined. 

_ 2*£M 

v Then, in step 614, cylinder exhaust gas recirculation amount is 
determined by filtering m egr . In particular, according to the 

present invention, the value is filtered synchronously with 
engine firing events. The filter coefficient (a) is a function 
20. of engine speed. The following equation shows the filtering' 
method . 

Using the embodiment of Figure 3, this can be rewritten as: 
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ik4map*dp . , 

m ^ = a T~N + (1 " a)m ^ 



l cyl L 

Next, in step 616, cylinder air amount is determined by 
subtracting the cylinder exhaust gas recirculation amount from 
the total cylinder amount as shown below. 

5 m cyt air =m cyt- m cyt egr 

Then, this value is used to calculate (open-loop) fuel injection 
amount (f) based on desired air-fuel ratio (a/f d ) . Also, fuel 
injection amount (f) can be adjusted based on a measured exhaust 
gas oxygen concentration from a HEGO sensor using methods known 
10 to those skilled in the art to provide closed loop air-fuel ratio 
control . 



alf d 



Thus, according to the present invention, it is possible to 
use the improved EGR flow estimation provided by having a 

15 downstream orifice and an upstream valve wherein differential 
pressure across the downstream orifice and manifold pressure 
provide the EGR flow estimate. Further, this improved EGR flow 
estimate is then filtered to account for manifold dynamics and 
used to calculate an improved cylinder air amount. Then, this 

20 improved cylinder air amount is used in air-fuel ratio control. 
In another embodiment, internal combustion engine 10 
comprising a plurality of cylinders, one cylinder of which is 
shown in Figure 7, is controlled by electronic engine controller 
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12. Exhaust manifold 48 is shown coupled to exhaust gas 
recirculation valve 70 via exhaust gas recirculation tube 202.. 
Exhaust gas recirculation valve 70 is also coupled to intake 
manifold 44 via orifice tube 74. Orifice tube 74 has orifice 205 
5 for restricting flow therein. Alternatively, the engine may be 
configured such that the fuel is injected directly into the 
cylinder- of the engine, which is known to those skilled in the 
art as a direct injection engine. Two-state exhaust gas oxygen 
sensor 98 is shown coupled to exhaust manifold 4 8 downstream of 

10 catalytic converter 20. Sensor 16 provides signal EGO to 

controller 12 which converts signal EGO into two-state signal 
EGOIS. A high voltage state of signal EG01S indicates exhaust 
gases are rich of a reference air/fuel ratio ctnd a low voltage 
state of converted signal EGOl indicates exhaust gases are lean 

15 of the reference air/fuel ratio. Sensor 98 provides signal EG02 
to controller 12 which converts signal EG02 into two-state signal 
EG02S. A high voltage state of signal EG02S indicates exhaust 
gases are rich of a reference air/fuel ratio and a low voltage 
state of converted signal EG02S indicates exhaust gases are lean 

20 of the reference air/fuel ratio. 

Controller 12 is shown receiving various signals from 
sensors coupled to engine 10, in addition to those signals 
previously discussed, including: a measurement of manifold 
pressure (MAP) from manifold pressure sensor 116 coupled to 

25 intake manifold 44; and a measurement of exhaust gas 
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recirculation pressure (EGRP) from exhaust pressure sensor 119 
coupled to orifice tube 74 upstream of orifice 205. In a 
preferred aspect of the present invention, engine speed sensor 
119 produces a predetermined number of equally spaced pulses 
5 every revolution of the crankshaft. 

Referring now to Figure 8, a flowchart of a routine 
performed by controller 12 to generate fuel trim signal FT is now 
described. A determination is first made whether closed-loop 
air/fuel control is to be commenced (step 122) by monitoring 

10 .engine operation conditions such as temperature. When closed- 
loop control commences, signal EG02S is read from sensor 98 (step 
124) and subsequently processed in a proportional plus integral 
controller as described below. 

Referring first to step 126, signal EG02S is multiplied by 

15 gain constant GI and the resulting product added to products 
previously accumulated (GI * EG02Si_i) in step 128. Stated 
another way, signal EG02S is integrated each sample period (i) in 
steps determined by gain constant GI . During step 132, signal 
EG02S is also multiplied by proportional gain GP . The integral 

2 0 value from step 12 8 is added to the proportional value from step 
132 during addition step 134 to generate fuel trim signal FT. 

The routine executed by controller 12 to generate the 
desired quantity of liquid fuel delivered to engine 10 and 
trimming this desired fuel quantity by a feedback variable 

25 related both to sensor 98 and fuel trim signal FT is now 
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described with reference to Figure 9. During step 158, an open- 
loop fuel quantity is first determined by dividing the difference 
between inducted mass air flow (AMPEM, created from the signal 
FMAP and RPM as described later herein with particular reference 
5 to Figure 10) , which includes both fresh charge and exhaust gas 
recirculation, and exhaust gas recirculation estimate (EM), which 
■is described later herein with particular reference to Figure 10, 
by desired air/fuel ratio AFd which is typically the 
stoichiometric value for gasoline combustion. However, setting 

10 AFd to a rich value will result in operating the engine in a rich 
state. Similarly, setting AFd to a lean value will result in 
operating the engine in a lean state. Also, signal AMPEM is 
constructed from FMAP and RPM in the common speed density method 
known to those skilled in the art and can be easily empirically 

15 determined. This open- loop fuel quantity is then adjusted, in 
this example divided, by feedback variable FV. 

After determination that closed- loop control is desired 
(step 160) by monitoring engine operating conditions such as 
temperature' (ECT) , signal EG01S is read during step 162 . During 

20 step 166, fuel trim signal FT is transferred from the routine 
previously described with reference to Figure 8 and added to 
signal EG01S to generate trim signal TS . 

During steps 170, 172, 176, and 178, a proportional plus 
integral feedback routine is executed with trimmed signal TS as 

25 the input . Trim signal TS is first multiplied by integral gain 
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value KI (step 170) , and the resulting product added to the 
previously accumulated products (step 172) . That is, trim signal 
TS is integrated in steps determined by gain constant KI each 
sample period (i) during step 172. A product of proportional 
5 gain KP times trimmed signal TS (step 176) is then added to the 
integration of KI * TS during step 178 to generate feedback 
variable FV. 

Calculating exhaust gas recirculation estimate (EM) is now 
described with particular reference to the diagram shown in 

10 Figure 10. In particular, Figure 10 shows how the upstream 
pressure (pi) , which is signal EGRP in this example, and 
downstream pressure (p2) , signal MAP in this example, are 
processed to form the signal EM. First, in block 1000, upstream 
pressure pi is processed through a first filter known to those 

15 skilled in the art as an anti-aliasing filter with a cut-off 
frequency equal to fl. Similarly, in block 1002, downstream 
pressure p2 is processed through a second anti-aliasing filter 
with a cut-off frequency equal to f 2 . In some applications, it 
is unnecessary to use either the first or the second anti- 

20 aliasing filter because the geometry of the exhaust gas 

recirculation creates a mechanical filter that removes the 
unwanted high frequencies. Further frequencies fl and f2 are set 
considerably higher than the necessary control bandwidth. 
Next, in block 1004, the result of block 1000 is 

25 synchronously sampled with an engine rotation signal, such as, 
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for example, RPM, such that the sampling is at a rate 
proportional to the firing frequency of the engine. For example 
the sampling rate could be twice the firing frequency of the 
engine. The proportion is generally chosen such that the 
5 sampling is at a rate of twice the highest harmonic frequency 
that contains significant energy. Also, as would be obvious to 
one of ordinary skill in the art and suggested by this 
disclosure, any multiple of firing frequency greater than that 
determined above could be used. If, for example, the exhaust gas 

10 recirculation and engine geometry are such that higher order 

harmonics are present in the upstream pressure signal pi, such 
as, for example, harmonics of twice or four times the firing 
frequency, a sampling rate of four or eight times the firing 
frequency may be necessary. Similarly, in block 1006, the result 

15 of block 1002 is synchronously sampled with engine speed signal 
RPM, such that the sampling is at a rate proportional to the 
firing frequency of the engine. Additionally, it is not 
necessary that the sampling rate be equal in blocks 1004 and 
1006. For example, block 1004 could synchronously sample at 

20 twice the firing frequency of the engine and block 1006 could 
sample at eight times the firing frequency of the engine. 

Alternatively, as is obvious to one of ordinary skill in the 
art and suggested by this disclosure, the pressure signal could 
be sampled at a frequency substantially proportional to the 

25 dominant frequency contained in the signal. This dominant 
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frequency is usually equal to firing frequency. Thus, sampling 
. at a rate proportional this dominant frequency could be 
accomplished using a circuit known to those skilled in the art as 
a phase-locked loop. However, because the phase locked loop 
5 shcme is sometimes searching for the dominant frequency during 
transients, this process may be suspended based on a change of 
position in throttle plate 62. During the transition, an open 
loop estimate of how the change in throttle plate 62 affects 
exhaust gas recirculation and manifold pressure must be obtained. 

.10 This can be done using a predetermined map obtained through 

testing or analytical procedures and is known to those skilled in 
the art, where the transient behavior is estimated based on 
change of position in throttle plate 62 and other operating 
conditions, such as for example engine speed. 

15 Next, digital filters in blocks 1008 and 1010 process the 

results of blocks 1004 and 1006. The digital filters, 
represented by G(z) or G' (z) used in blocks 1008 and 1010 are 
known to those skilled in the art as digital notch filters. In 
this application, each notch filter removes the firing frequency 

2 0 (and higher harmonics if necessary) of the engine. The equation 
below represents an example of a notch filter in the discreet 
domain for sampling at a rate of twice the firing frequency Use 
of a notch filter G(z) is also described later herein with 
particular reference to Figure 5. 

25 
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G(z) = (1+ z" 1 ) /2 
If the sampling were done at a rate of eight times the firing 
frequency, then the following notch filter would be used as 
described by G ' ( z) . Again, while this removes unwanted 
5 frequencies, transient performance is not hindered. Use of a 
notch filter such as G' (z) is described later herein with 
particular reference to Figure 6 and 7. 

G 1 (z) = (1+ z" 1 + z' 2 + z~ 3 + z" 4 + z~ 5 + z" 6 + z" 7 )/8 

10 

The digital filter may be different between blocks. 1008 and 1010 
and different than that shown above if necessary, such as if, for 

example, the geometry of the exhaust gas recirculation system was 

-i - 
such that the certain frequencies were excessively amplified due 

15 to resonances. Also, the filter may be different between blocks 
1008 and 1010 if block 1004 synchronously sampled at twice the 
firing frequency of the engine and block 1006 sampled at eight 
times the firing frequency of the engine. 

The pressure difference is then created by subtracting the 

20 output of block 1010, which is filtered manifold pressure FMAP, 
from the output from block 1008 . This pressure difference is 
then used in block 1012 to create signal EM through a 
predetermined map or equation between pressure difference and 
exhaust gas recirculation flow, and, if necessary, engine 

25 operating conditions. For example, exhaust gas temperature may 
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be used to adjust the calculation of exhaust gas recirculation 
flow. 

Also, in block 1014, signals FMAP and RPM are used to 
calculate the mass of gas flow entering the cylinder (AMPEM) . 
5 The common speed density equations known to those skilled in the 
art are used to convert the filtered manifold absolute pressure 
with the engine speed to the total mass of gas (exhaust gas and 
fresh air charge) entering the cylinder. If necessary, these 
basic equations can be modified by engine operating conditions, 

10 such as for example gas temperature, or any other condition known 
to those skilled in the art and suggested by this disclosure. 

Thus, an estimate of the exhaust gas recirculation and fresh 
air entering the cylinder is obtained that is substantially free 
of unwanted frequencies yet retains a bandwidth that is much 

15 greater than would be obtained with conventional filtering 

methods. Thus, the estimate can more accurately track transient 
operation and yield more accurate air/fuel ratio control. 

An example of synchronously sampling a waveform is now 
described with particular reference to the plot shown in Figure 

2 0 11. A fluctuating pressure signal, shown by the solid line and 
labeled A, is sampled with a frequency equal to twice the 
frequency of the actual signal. The sampled values are shown by 
points. The reconstructed waveform based on the synchronously 
sampled values and the filter previously described 'herein with 

25 particular reference to the function G(z) is shown as the dotted 
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line and labeled B. For comparison, a signal using a 
conventional low pass filter, which is required for conventional 
sampling schemes, is shown by a dash dot line and labeled C. In 
this example, the exhaust gas recirculation estimate formed using 
5 the synchronous sampling will yield a more accurate value that 
will allow for better overall air/fuel ratio control. 

Another example of synchronously sampling a waveform is now 
described with particular reference to the plot shown in Figure 
12. A fluctuating pressure, shown by the solid line and labeled 

10 D, is sampled with a frequency equal to eight times the frequency 
of the lowest harmonic order. This signal represents a typical 
exhaust pressure during steady state operating conditions. The 
sampled values are shown by points. The reconstructed waveform 
based on the synchronously sampled values and the filter 

15 previously described herein with particular reference to the 

function G' (z) is shown as the dotted line and labeled E. This 
result could not be obtained unless the sampled values are all 
perfectly spaced with the rotation of the engine, the synchronous 
sampling frequency was such that it was twice the highest 

2 0 significant harmonic frequency of the pressure signal, and the 

appropriate notch filter was used. In this example, the air flow- 
entering the cylinder estimate formed using the synchronous 
sampling will yield an accurate value that will allow for optimal 
air/fuel ratio control. 
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Now referring to Figure 13A-13B and in particular to Figure 
13A, the plot shows the frequency content of the pressure 
waveform shown in Figure 12. This pressure could represent, for 
example, the exhaust manifold pressure for a steady state firing 
5 frequency of the engine of approximately 50 Hz. Figure 13B shows 
a plot of the magnitude versus frequency of the filter G' (z) . 
Thus, the scheme previously described herein with particular 
reference to Figure 10, comprises (in the frequency domain) 
multiplying the plots of Figures 13A and 13B. This shows that 
...10 the mean value, or DC component as known to those skilled in the 
art, is preserved. The result is a signal substantially free of 
undesirable frequencies for mean value model computations. 
There are also other alternative embodiments of the present 
.invention. For example, using a synchronous sampling scheme is 

15 not dependent on the orifice being located downstream of the 

exhaust gas recirculation flow control valve. The scheme could 
be employed using a pressure sensor upstream and a pressure 
sensor downstream of the orifice, with the exhaust gas 
recirculation flow control valve still between the downstream . 

20 pressure sensor and the intake manifold, as in current production 
vehicles. Furthermore, the method is not restricted to flow 
measurement with an orifice. Other flow measurement techniques 
known to those skilled in the art could be used with the above 
described method such as, for example, a venture, a pi tot tube, 

25 or a laminar flow element. 
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This concludes the description of the Preferred Embodiment. 
The reading of it by those skilled in the art would bring to mind 
many alterations and modifications without departing from the 
spirit and scope of the invention. For example, rather than a 
5 differential pressure sensor and an absolute downstream pressure 
sensor, two absolute pressure sensors (one upstream of the fixed 
area orifice ( p x ) and one downstream of the fixed area orifice 
■A Pi)) could be used as shown above. From the two absolute 
pressure sensors, a differential pressure signal and a downstream 
10 pressure signal could be produced. In other words, manifold 

pressure could be provided from the downstream pressure sensor 
and the differential pressure could be provided from the 
difference between the upstream and downstream pressure sensors. 
-The difference could be calculated internally in the 
15 microprocessor or in an analog circuit coupled to the two 

pressure sensors. In other words, the equation of step 216 would 
become : 

EM = Kj( Pl -p 2 )*p 2 

Accordingly, it is intended that the scope of the invention be 
20 defined by the following claims. 
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